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COMPUTER STUDY OF A DEAD-BAND CONDENSING-PRESSURE 
CONTROL FOR SNAP-8 STARTUP 
by Roy C. Tew and James D. Gallagher 
Lewis Research Center 

SUMMARY 

A digital-computer simulation of the SNAP-8 space power system was used to eval- 
uate a dead -band method of controlling condensing pressure during the startup transients. 
With this control method, the mercury-condensing pressure is controlled by varying the 
condenser -coolant flow rate whenever pressure goes outside some predetermined dead- 
band. Increasing the coolant flow, for example, tends to increase the rate of heat ex- 
traction from the condenser and therefore to decrease the condensing pressure. Since 
the condensing pressure is also the turbine back -pressure, high pressures must be pre- 
vented so as not to excessively lower turbine output power. And, since the mercury 
pump net positive suction head (NPSH) is maintained by the condensing pressure and the 
amount of subcooling, it is equally important to avoid excessively low pressure through- 
out startup in order to prevent the pump from cavitating. 

The study covered a range of values for each of the following control parameters: 

(1) pressure dead-band width, (2) rate of change of coolant flow, and (3) initial coolant 
flow through the condenser. The study also covered a range of values for each of the 
following system characteristics: (1) the amount of mercury that accumulates in the 
condenser (condenser inventory); and (2) the mercury -flow ramp rate, which determines 
the rate of buildup of the mercury liquid -vapor interface in the condenser. 

In addition to detailed information on the control characteristics, the study produced 
the following results concerning the sensitivity of control performance to system char- 
acteristics: 

(1) There is a limited range of inventories for which control performance is satis- 
factory. 

(2) Mercury -flow ramps to the system self-sustaining level from 40 to 100 seconds 
in duration should allow satisfactory control performance (provided condenser inventory 
is within the correct range). 



INTRODUCTION 


SNAP-8 is a 35-kilowatt electric nuclear Rankine-cycle space power system cur- 
rently under development. It has three liquid-metal heat -transfer loops and a lubricant- 
coolant loop. (A schematic diagram showing the three liquid metal heat transfer loops 
is shown in fig. 1.) In the first loop NaK, the eutectic mixture of sodium and potassium, 
is circulated to transfer heat from the nuclear reactor to the mercury boiler. The 
second loop is the mercury Rankine loop. In this loop, liquid mercury is heated, evap- 
orated, and superheated in a counterflow multitube boiler. After passing through the 
turboalternator, the mercury is condensed and subcooled in a counterflow multitube 
condenser. The third loop, the heat-rejection loop, circulates NaK to transfer heat 
from the condenser to the radiator, where it is dissipated to space. The fourth loop 
circulates oil in order to cool and lubricate the turboalternator and the motor-driven 
mercury pump. This electric generating system is being developed to operate unattended 
for at least 10 000 hours in space after automatic startup. It will also have shutdown 
and restart capabilities. 

The startup of SNAP-8 has two phases: (1) reactor startup during Which the primary 
loop is brought to design temperature; and (2) power -conversion system startup during 
which self-sustaining operation and, subsequently, full-power operation is attained. At 
the beginning of the second phase, the mercury inventory is injected into the evacuated 
second loop in accordance with a programmed ramp of flow rate as a function of time so 
that the turboalternator is brought to rated speed. When the injection process is com- 
plete, the mercury pump recirculates the liquid mercury from the condenser back to the 
boiler. When the system reaches the self-sustaining power level, the mercury flow rate 
is held constant until transients have settled out; it is then increased in a gradual manner 
to the full-power value. 

During power conversion system startup, the condensing pressure should vary 
smoothly from essentially zero to the design value. There may be transient conditions, 
however, which will tend to cause the pressure to rise much higher than the design 
value (14. 0 psi, 9. 66 N/cm ). These high pressures will occur if the NaK flowing 
through the condenser is too low relative to the heat input from the mercury loop. And, 
since the condensing pressure is also the turbine backpressure, pressures much higher 
than 15. 0 or 16. 0 psi (10. 3 or 11. 0 N/cm ) should be avoided so as not to lower the tur- 
bine output power, especially when the turbine is accelerating. A very high condensing 

p 

pressure could also exceed the condenser structural design limit (60. 0 psi, 41. 4 N/cm ). 
It is also important to have sufficient condensing pressure throughout startup to prevent 
the mercury pump from cavitating. This is especially important under zero-gravity con- 
ditions, for which the pump net positive suction head (NPSH) is maintained by a com- 
bination of condensing pressure and subcooling. Too low a pressure will reduce the 
pump suction head below the level required to prevent cavitation. Therefore, some form 
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(a) Initial conditions for mercury loop startup. 
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(b) Conditions at design operating point. 
Figure 1. - System conditions. 
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of pressure control is required to avoid the operating problems associated with these 
two extremes of condenser pressure. 

This study is an evaluation of a dead -band method of controlling the mercury- 
condensing pressure during startup. With this method, the condensing pressure is con- 
trolled by varying the NaK coolant flow rate through the condenser whenever the con- 
densing pressure goes outside some predetermined dead-band. If, for example, the 
condensing pressure exceeds the upper dead-band limit, the coolant flow increases. 

The increased coolant flow increases the rate of heat extraction from the condenser 
and tends to reduce the condensing pressure. The coolant flow continues to increase 
until the pressure drops below the upper dead-band limit. 

In reference 1, condensing-pressure control for long-term operation of the SNAP-8 
system was studied. This study indicated that variations of condensing pressure around 
the design operating point could be expected, and methods of controlling these variations 
were investigated. One of the methods investigated was a dead-band control, and the 
results indicated that this method of control could be used successfully. These studies, 
however, were made only for pressure variations around the design operating point. 

The analysis herein treats the more critical variations occurring during the startup 
transient. 

The study was made by using a digital-computer simulation of the SNAP -8 system; 
it covered a range of values for each of the following control parameters: (1) pressure 
dead-band width, (2) a valve gain defined as rate of change of flow over flow, and (3) 
initial coolant flow through the condenser. Since the inventory in the condenser strongly 
affects the condensing pressure, runs were made for a range of condenser inventories 
(i.e. , by varying the amount of mercury injected into the system). The effect of mer- 
cury flow rate schedule was investigated by using two different mercury ramps to the 
self-sustaining level. 


DESCRIPTION OF DIGITAL COMPUTER SIMULATION 

The following major components of the SNAP-8 system were included in the digital 
computer dynamic simulation: boiler, turbine, alternator (and the associated electrical 
components), condenser, radiator, auxiliary start heat exchanger, heat-rejection loop, 
auxiliary start loop, flow-control valve, and all the pipe lines between these components. 
A reactor simulation was not used; therefore, the temperature of the NaK coming into 
the boiler was assumed to stay at the design value for all startup runs evaluated. (Ex- 
pected variations in boiler inlet temperature during startup should not significantly 
affect condenser input variables . ) 


4 



Condenser Simulation 


In the SNAP-8 condenser, the mercury flows axially in 73 inner tubes, while the 
NaK flows through the shell counterflow to the mercury. The tubes are tapered to 
maintain vapor velocity which provides a continual movement of condensate to the liquid- 
vapor interface in zero-gravity operation. A detailed description of this condenser and 
its steady-state performance is given in reference 2. Since the results used in this 
study depend heavily upon the condenser simulation, this simulation is described in 
detail in appendix A. 


Description of Control Logic 

Prior to mercury injection there is an initial NaK flow rate through the condenser, 
and the flow through the auxiliary start loop is at a value necessary to remove heat from 
the first loop for reactor startup. The auxiliary start loop includes the NaK flow path 
from the flow control valve to the auxiliary heat exchanger and back to the heat rejection 
loop (see fig. 1). During the preinjection period NaK flows are maintained by running 
the pumps on a battery-powered inverter at half design speed. The pumps are brought 
to design speed by switching them from the inverter to the alternator when the turbine 
speed passes through its half design speed level; thus, the pumps are said to be boot- 
strapped to their design level. In the simulation, the primary loop, heat-rejection 
loop (HRL), and auxiliary start loop (ASL) flows are bootstrapped at a rate proportional 
to the turbine speed from the time the turbine passes through its half design speed level. 
Therefore, the initial NaK flow at design pump speed is two times the initial NaK flow 
before bootstrapping. The flows are then kept at these levels during mercury injection 
until the condensing pressure exceeds some predetermined upper limit (the upper dead- 
band limit). At this time, both HRL and ASL flows are varied simultaneously in a 
manner such that the percentage change in flows with respect to time is constant. The 
flow through the condenser is increased as the ASL flow is decreased until either the 
condensing pressure drops below the upper dead-band limit or the design flow rate is 
reached. The flows do not change when the condensing pressure is within the dead-band. 
However, if the pressure goes below some predetermined lower limit (the lower dead- 
band limit), the flow through the condenser is decreased and the ASL flow is increased 
until either the pressure goes back above the lower dead-band limit or the condenser 
flow reaches the initial value for design pump speed. Therefore, this initial value of 
condenser flow is also the minimum flow; it is determined by the setting of the flow- 
control valve prior to mercury injection. (The flow returned to the initial (minimum) 
value for only one of the runs discussed. ) 
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The simulated flow -control valve produces the same rate of change of flow over flow 
regardless of the flow rate; that is, the ratio (1/w) x (dw/dt) is constant, but the differ- 
ential dw/dt (rate of change of flow rate) increases as the flow increases. The total 
heat -rejection-loop flow is the sum of the condenser coolant flow and the ASL flow. The 
rate of change of flow over flow for the condenser coolant flow and for the ASL flow are 
equal in magnitude but opposite in sign. The magnitude of this flow parameter can be 
considered the flow -control -valve gain; it will be referred to as such from now on. 

The value of the condensing pressure is obtained from the condenser simulation and 
is used together with a first-order lag (with a time constant of 0. 1 sec, which was in- 
cluded to simulate pressure sensor dynamics). 


Procedures for Making Startup Runs 

The initial conditions for all three loops, prior to mercury loop startup, are shown 
in the schematic diagram in figure 1(a) for condenser NaK flow equal to 9 percent of the 
design flow. For the case shown, the initial condenser flow (and the minimum flow) after 
bootstrapping of the pumps to design speed would be 18 percent of the design flow. In the 
first loop the NaK flow rate is at 50 percent of design value, the boiler is at a uniform 
temperature of 1300° F, and the auxiliary -start heat exchanger is extracting heat from 
the primary NaK loop. In the third loop the temperatures are determined by the NaK 
flow rate, the heat coming from the auxiliary start loop, and the radiator characteris- 
tics. The flow through the condenser keeps it at a uniform temperature equal to the ra- 
diator exit temperature . System conditions at the design operating point are shown in 
figure 1(b). For one group of startup runs, the mercury flow rate schedule (as seen at 
the boiler inlet) was as follows: (1) a ramp from zero to 53. 5 percent of design flow rate 
in 100 seconds; (2) remained at 53. 5 percent of design flow. This flow rate schedule is 
shown in figure 2(a), and it will be referred to as the slow injection schedule or the slow 
ramp rate from now on. It should be pointed out that 53. 5 percent of design flow rate is 
the minimum value needed for self-sustaining operation of the system; that is, operation 
with enough turboalternator electrical power to run the pumps and control components. 
For other startup runs the mercury was ramped from zero to 53. 5 percent of design 
flow rate in 40 seconds, and then remained at 53. 5 percent until die end of the run. This 
injection schedule will be referred to as the fast injection schedule or the fast ramp rate 
from now on (fig. 2(b) solid curve). The dashed curve in figure 2(b) indicates the form 
(but not necessarily the correct time intervals) that the complete mercury flow schedule 
is expected to have during startup of the SNAP-8 system. The slow ramp from 53. 5 
percent of design flow to design flow was not covered in the computer study. The slow 
transients resulting from this ramp were not expected to present a problem for the 
control and would, if studied, require a relatively large amount of computer time. 
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(b) Fast injection schedule. 

Figure 2. - Mercury injection schedules for startup. 

For a given mercury flow rate schedule, the injection process continued until the 
mercury injected was sufficient to give the desired condenser inventory. At the same 
time that the injection process was stopped, the isolation valve (fig. 1) was opened to 
allow mercury to flow from the condenser to the mercury pump. The mercury injection- 
isolation valve was located at the exit to the condenser. The line between this valve and 
the pump, the pump, and the line between the pump and the flow-control valve (which was 
located at the boiler inlet) were assumed to have been prefilled. 

RESULTS AND DISCUSSION 
Effect of Pressure Dead-Band Width and Valve Gain 

Pressure dead-band width and valve gain are considered together because both var- 
iables significantly affect the stability of the control. It is desirable to have the smallest 
dead-band and the highest gain consistent with stable operation of the control. Decreasing 
the dead -band width increases confidence that the condensing pressure is within an ac- 
ceptable range (within the limits of stable operation). However, decreasing the dead- 
band width also makes the response of condensing pressure more oscillatory. Increasing 
the valve gain makes the control respond faster (thus decreasing overshoot of the dead- 
band limits) but also makes condensing-pressure response more oscillatory. The runs 
discussed allow approximate determination of the smallest dead-band width and the high- 
est valve gain consistent with stable operation of the control. 
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(a) Pressure dead-band width, 2 psi (1.4 N/cm 2 ). 



(b) Pressure dead-band width, 3 psi (2. 1 N/cm 2 ). 



Time, sec 

(c) Pressure dead-band width, 4 psi (2.8 N/cm 2 ). 

Figure 3. - Effect of dead-band width on control response. Valve gain, 
3. 2 percent per second; initial condenser coolant flow, 18 percent 
of design flow (after bootstrapping); final condenser inventory, 

40 pounds mass (18.2 kg); slow mercury flow schedule (100-sec 
ramp). 
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For all startup runs discussed in this section, the slow mercury flow schedule (100- 
sec ramp) was used and the final condenser inventory was 40 pounds (18 kg). 

In figure 3 mercury-condensing-pressure variations for startup runs with three 
different pressure dead-hand widths are shown. For all three runs the valve gain was 

3. 2 percent change in flow per second and the NaK flow, immediately after bootstrapping, 
was 18 percent of the design flow. (The NaK flow was 9 percent of design flow before 
bootstrapping. ) Initial conditions for the three runs were the same. Figures 3(a) to (c) 
correspond to dead-band widths of 2.0, 3.0, and 4.0 psi, respectively. A trend toward 
less oscillatory response is seen in going from the 2.0-psi (1. 4-N/cm ) dead-band run 

9 

to the runs with larger dead-band widths. The 2. 0-psi (1. 4-N/cm ) dead-band run 

p 

shows a considerably more oscillatory response than the 3. 0-psi (2. 1-N/cm ) run. The 
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4. 0-psi (2. 8-N/cm ) dead-band run was slightly less oscillatory than the 3. 0-psi (2. 1- 
N/cm ) run; all remaining runs were made using a 4. 0-psi (2. 8-N/cm ) dead-band. 

In figure 4 mercury-condensing-pressure plots for startup runs with three different 
NaK flow-control valve gains are shown; the pressure dead-band width used was 4. 0 psi 
(2. 8 N/cm ). All other parameters and the initial conditions are the same as for the 
runs of figure 3. Figure 4(a) to (c) correspond to valve gains of 3. 2, 6.4, and 12-. 8 per- 
cent per second, respectively. A trend toward more oscillatory response of condensing 
pressure is seen in going from the lowest-gain run to the runs with larger gains. The 
run made with a valve gain of 12.8 shows a much more oscillatory response than the 
6. 4-gain run. The 3. 2 run is slightly less oscillatory than the 6. 4 run, but shows more 
overshoot of the upper dead-band. It might be desirable to tolerate more oscillations 
with the 6.4 gain in order to get less overshoot than that observed with the 3.2 gain. The 
remaining runs discussed in this report were made using a valve gain of 6. 4 percent per 
second. 

The six plots of condensing pressure previously discussed are identical up to the 
point at which the pressure exceeds the upper dead-band limit for the first time. 


Effect of Initial NaK Flow 

The control parameter, initial condenser NaK flow, could be made high enough to 
eliminate the initial overshoot. If this parameter is made too high, however, the con- 
densing pressure (and the mercury pump NPSH) will be too low at the end of mercury 
injection. The runs discussed allow an approximate determination of the highest initial 
condenser NaK flow consistent with adequate NPSH at the end of mercury injection (for 
the given Hg ramp rate and condenser inventory). 

In figure 5 mercury-condensing pressure for startup runs with three different ini- 
tial NaK flows are shown; the valve gain used was 6. 4 percent change in flow per second. 
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(a) Valve gain, 3.2 percent per second. 



(b) Valve gain, 6.4 percent per second. 



Time, sec 

(c) Valve gain, 12.8 percent per second. 

Figure 4. - Effect of valve gain on control response. Pressure 
dead-band width, 4 psi (2.8 N/cm^) ; initial condenser coolant 
flow, 18 percent of design flow (after bootstrapping); final 
condenser inventory, 40 pounds mass (18.2 kg); slow mer- 
cury flow schedule (100-sec ramp). 
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Mercury-condensing pressure, N/cm 



(a) Initial condenser coolant flow, 18 percent of design flow (after bootstrapping). 



(b) Initial condenser coolant flow, 20 percent of design flow (after bootstrapping). 



0 50 100 150 200 250 300 350 


Time, sec 

(c) Initial condenser coolant flow, 25 percent of design flow (after bootstrapping). 

Figure 5. - Effect of initial condenser coolant flow on control response. Valve 
gain, 6.4 percent per second; pressure dead-band width, 4 psl (2.8 N/cm 2 ); 
final condenser inventory, 40 pounds mass (18.2 kg); slow mercury flow 
schedule (100-sec ramp). 
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All other parameters are the same as for figure 4. The heat-rejection-loop initial 

temperatures were not changed when the initial NaK flow was changed. The initial NaK 

flow mentioned here corresponds to the initial HRL temperatures and is one -half of the 

"initial” NaK flow at design pump speed. (In a real situation, maintaining the same 

radiator exit temperature and increasing the flow through the condenser from 9 to 

10 percent of design while maintaining a 38 -kW heat-rejection level implies a drop in 

radiator inlet temperature from 230° to 222° F (373 to 369 K). These conditions imply 

a radiator 2 percent larger in area.) Therefore, the initial conditions used, with the 

exception of the initial NaK flow, were the same as for the previous runs. Figures 5(a) 

to (c) correspond to initial NaK flows (at design pump speed) of 18, 20, and 25 percent of 

design HRL flow, respectively. Changing the initial NaK flow from 18 to 20 percent of 

the design flow reduces the pressure overshoot of the upper dead-band. However, the 

pressure level at the end of mercury injection is smaller for the 20 -percent than for 

the 18 -percent initial flow case. A further increase in initial NaK flow to 25 percent of 

the design value reduces the condensing pressure level at the end of injection to about 

2 

0. 78 psi (0. 54 N/cm ). From figure 6 the required minimum NPSH at 5.3. 5 percent of 

o 

design mercury flow is 1.0 psi (0.69 N/cm ). Therefore, a condensing pressure this 
low should be expected to result in mercury pump cavitation. The 20 -percent initial 

Q 

NaK flow case yields a condensing pressure of 2. 36 psi (1. 63 N/cm ) at the end of 
mercury injection. This pressure is 1. 36 psi (0.94 N/cm ) above the minimum required 
NPSH and is therefore adequate. However, to provide more NPSH margin at the end of 
injection it may be desirable to tolerate slightly more overshoot and use an initial NaK 
flow of 18 percent of design. For the 18 -percent initial NaK flow case, the pressure at 
the end of injection is 3.2 psi (2.2 N/cm ) above the minimum required NPSH. 



Figure 6. - Minimum required net positive suction head (NPSH) for 
mercury pump as function of design mercury flow rate. 
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Effect of Mercury Condenser Inventory and Mercury Flow Rate Schedule 


Condenser mercury inventory and the mercury flow rate schedule are system char- 
acteristics which are determined primarily by conditions outside the condenser. Con- 
denser inventory is uncertain, even when a fixed inventory is injected into the system, 
for the following reasons: 

(1) Condenser mercury inventory is dictated by boiler inventory. Boiler heat- 
transfer rate, and therefore boiler inventory, is very sensitive to contaminants on the 
inner surfaces of the boiler tubes. 

(2) When the system is restarted, the amount of mercury remaining in the system 
since the previous shutdown cannot be accurately determined. 

Mercury ramp rate is dictated primarily by other system considerations. Ramps of 
duration much less than 40 seconds do not appear advisable from the viewpoint of reac- 
tor temperature transients. Very long ramps (150 sec or longer) can result in the 
application of electrical load (pumps, control components) while mercury flow, and 
therefore applied turbine torque, is too low to accomplish bootstrapping. Even so, 
acceptable range of these inputs may be limited by condenser -pressure control require- 
ments. The runs discussed in this section allow evaluation of the control performance 
over ranges of mercury inventory and mercury flow ramp rate. 

Runs made with ranges of condenser inventories are shown in figures 7 and 8. The 
slow mercury flow rate schedule (100-sec ramp) was used for the runs in figure 7 and a 
fast mercury flow rate schedule (40 -sec ramp) was used for the runs shown in figure 8,. 
Initial NaK flow rate through the condenser of 20 percent of design flow (after boot- 
strapping) and a flow-control -valve gain of 6. 4 percent change in flow per second were 
used. 

The condensing pressure during startup is seen to be very sensitive to changes in 

condenser inventory. For the runs shown in figure 7 (100 -sec Hg ramp) the pressure 

at the end of injection is seen to be 0. 33 psia (0. 23 N/cm ) for the smallest inventory of 

17 pounds (7. 7 kg). Increasing the inventory to 40 pounds (18 kg) raises the pressure at 

o 

the end of injection to 2.4 psia (1. 6 N/cm ) as seen in figure 9(b). Condensing pressure 
at the end of injection is shown as a function of condenser inventory in figure 9 for both 
the 100-second ramps of figure 7 and for the 40-second ramps of figure 8. 

From figures 7 to 9 it is possible to determine acceptable ranges of condenser in- 
ventory for the fast and slow mercury flow rate schedules. Approximate minimum in- 
ventories for each of the two ramp rates may be determined from figure 9. For the 
40 -second mercury ramp the end of the injection process always occurred after the 
mercury flow had reached 53. 5 percent of the design flow. From figure 6 the minimum 
required NPSH at 53. 5 percent of design flow is 1. 0 psia (0. 69 N/cm 2 ). For the 
smallest inventories used with the 100-second ramp, the injection process ended before 
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Mercury condensing pressure, N/cm' 



(c) Final condenser inventory, 53 pounds mass (24.0 kg). (d) Final condenser inventory, 64 pounds mass 129. 0 kg). 

** Figure 7. - Effect of condenser inventory on control response for slow mercury flow schedule. Valve gain, 6.4 percent per second; pressure 

dead-band width, 4 psi (2.8 N/cttt); initial condenser coolant flow, 20 percent of design flow (after bootstrapping); slow mercury flow schedule 
{100-sec ramp). 

completion of the mercury ramp. However, the lowest flow level at which injection 
ended was 45 percent of design. The minimum mercury pump NPSH required at 45 per- 
cent of design flow is 0. 75 psia (0. 52 N/cm^), or very close to 1.0 psia (0. 69 N/cm\ 
as seen from figure 6. Therefore, assuming a minimum required NPSH of 1. 0 psia 
(0. 6 N/cm ) for both ramp rates, minimum required inventories for both rates may be 
determined from figure 9. Minimum inventory for the 40 -second ramp is 25 pounds 
(11 kg) and for the 100 -second ramp is 30 pounds (14 kg). 

Inventory maximums are dictated by initial pressure overshoots of the upper dead- 
band. By comparison of figures 8(d) and (e), 60 pounds (27 kg) appears a reasonable 
approximation of the maximum allowable inventory for the 40 -second ramp. Figure 7(d) 
shows that a 64-pound (29 -kg) inventory would be acceptable for the 100 -second ramp 
and even a few more pounds (kg) would probably be acceptable. 

Therefore, for a 40-second mercury ramp a 35 -pound (16 -kg) range of inventories 
from 25 pounds (11 kg) to 60 pounds (27 kg) is acceptable. For a 100 -second mercury 
ramp at least a 35 -pound (30 -kg) range from 30 pounds (14 kg) to 65 pounds (30 kg) 
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Mercury-condensing pressure, psi 


5H 
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(a) Final condenser inventory, 22 pounds mass (10. 0 kg}. (b) Final condenser inventory, 32. 5 pounds mass (14. 5 kg). 




0 50 100 150 200 250 300 

Time, sec 


(e) Final condenser inventory, 66 pounds mass (30.0 kG). 

Figure 8. - Effect of condensing inventory on control response for fast mercury flow schedule. Valve gain, 6.4 percent per second (2. 8 N/cm 2 }; 
pressure dead -band width, 4 psi (2. 8 N/cm2) ; initial condenser coolant flow, 10 percent of design flow (before bootstrapping); fast mercury flow 
schedule (40-sec ramp). 
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Figure 9. - Condensing pressure at end of mercury injection as function of condenser in- 
ventory. 


appears acceptable. It also follows that any mercury ramp between 40 and 100 seconds 
would probably be acceptable for inventories in the 30-pound (14-kg) range from 
30 pounds (14 kg) to 60 pounds (27 kg). 


SUMMARY OF RESULTS 

The results from a study of a dead-band condensing -pres sure control for SNAP-8 
startup can be summarized as follows: 

(1) Control response becomes more oscillatory as dead-band width is decreased. 
Response goes from a non-limit-cycling to a limit-cycling form for widths between 
3. 0 psi (4. 1 N/cm^) and 2.0 psi (1. 4 N/cm^). 

(2) Control response becomes more oscillatory as the valve gain is increased. 
Response goes from a non -limit -cycling to a limit -cycling form for valve gains between 
6. 4 and 12. 8 percent per second. 

(3) Increasing the initial NaK flow decreases the initial overshoot but also decreases 
the mercury pump NPSH during startup. Mercury pump NPSH goes from adequate to 
inadequate when initial NaK flow is increased from 20 to 25 percent of design. An 

18 -percent initial NaK flow gives a more comfortable NPSH margin than a 20 -percent 
initial NaK flow and still has an acceptable initial pressure overshoot of the upper dead- 
band limit. 
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(4) There is a limited range of condenser inventories for which control performance 
is satisfactory. Satisfactory performance was obtained over a 35 -pound (16 -kg) range 
for both mercury flow ramps. The satisfactory range extended from 25 pounds (11 kg) 
to 60 pounds (27 kg) for the 40-second ramp and from 30 pounds (14 kg) to 65 pounds 

(29 kg) for the 100-second ramp. 

(5) It follows from (4) that mercury flow ramps to the system self-sustaining level 
from 40 to 100 seconds in length should yield satisfactory performance provided con- 
denser inventory is in the corresponding 35 -pound (16-kg) range. It also follows that if 
the inventory is restricted to the 30 -pound (14 -kg) range from 30 pounds (14 kg) to 

60 pounds (27 kg), all ramps from 40 to 100 seconds in length should yield satisfactory 
performance. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 20, 1969, 

120-27. 
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APPENDIX A 


SNAP-8 CONDENSER SIMULATION 

The SNAP-8 condenser is illustrated in figure 10; the simulation follows. Symbols 
are given in appendix B. For simplification the following assumptions are made: 

(1) Condenser heat transfer along the flow axis is negligible. 

(2) A mean diameter can be used for all calculations, which equals 1/2 (Diarm n | e< . + 
Dlan Wlet>' 

(3) The heat capacity of the tubes can be added to the heat capacity of the NaK. 

(4) Condensing takes place right from the beginning of the condenser (i.e. , no 
superheat region) . 

(5) The heat-transfer area is 73 times the heat-transfer area for one tube. 

With these assumptions, the condenser can be considered to be a simpler tube -in-tube 
condenser, as shown in figure 11. 

The condenser is divided into 40 equal segments or lumps (41 nodes) and is con- 
sidered to be made up of two main regions: the condensing region and the liquid region 




NaK in NaK out 

Figure 11. - Simplified condenser model. 
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(fig. 11). The average NaK temperature for the mercury-condensing region can be cal- 
culated as follows: 


where 


K=I 


NaK, av 


y T NaK,K 

A— — ■ 

Cond Ki 

K=1 

region 


Cond 

region 


Number of nodes 



(see fig. 12). This assumes that NaK temperature is constant from the last lump 


Interface location 



Figure 12. - Location of condenser 
liquid-vapor interface. 


before the interface to the interface. And, with assumption (4), T NaK , and the heat 
flow into the condenser from the turbine subroutine, Q- n , the condensing temperature 
can now be obtained. 


Qin= W T 


cond 


T NaK, a^ 


T 

cond 


Q 


in 


U c A c 


T NaK, av 
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It should be noted that A c = A c j n ^ a i x V c /V c ^ n ^j a V Therefore, as the condensing 
volume V c changes so does the condensing area A Q . 

Having obtained the condensing temperature, the condensing pressure can now be 
calculated from the Clausius -Clapeyron equation 


d(ln P) = _Ah_ 

d T rt 2 


T, P=constant 


where 

Ah heat of vaporization 
R specific gas constant 

This equation can be used because, for small temperature changes, Ah can be con- 
sidered constant. With this assumption, the preceding equation can be integrated (ref. 3) 
as follows: 



where 


C 


1 “ 


Ah 

RTq 


+ lnP Q 


C 


2 


Ah 

R 
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And for this case the condensing pressure becomes 


P cond ~ 


12 890 


T cond +460 


where P Q , T q correspond to a point on the mercury vapor saturation curve. 

The change in the condensing volume with time is calculated from the mass balance 
equation 


dV 


c _ 


AWtt„iv 


"Hg'Hg 


dV 

i 

dt 


c _ 


(w H g> ' J 


vapor 


Hg' 


liquidj 


Hg 


Approximating this differential over a finite time interval At yields 


d V 

(V ) = (V ) + — £ (At) 

L new c old d t 

^Hg( At ) 

The temperature at each node in the condensing region can be set equal to the condensing 
temperature 


(VJ = (VJ 

c new c old 


(W Hg ) 


vapor 


-(W H g) 


liquid 


T Hg, K “ T cond K = 

The temperature of the mercury at each node K 
as follows: 

8 [ p Hg (Cp) Hg T Hg,K] _ 9 [pHg V Hg (Cp) Hg T H g) K] , U Z (AT K )(dA HT ) 

at a x (dv £ ) 


1, 2, 3, . . I 

in the liquid region can be calculated 


Assuming p Hg , (Cp) H g, and v Hg are constant with respect to distance along the con- 
denser, and approximating this differential over a finite interval AX, which is equal to 


one lump, results in 
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dT Hg,K _ p Hg V Hg( Cp W( T Hg,K+l ~ T Hg, 
dt p Hg (Cp) Hg (AX) 


2U Z 7rr HT^ AX ^ AT AK^ 

{c P ) H ^y(AX) 


where 


K = 1+ 1, 1+ 2, • • • , 40 

AT AK = \ [f T NaK, K+l + T NaK,K^ " ^ T Hg,K+l " T Hg, K^] 

For numerical stability, AX should always be less than v Hg (At). When this is true, 
this differential equation can be approximated over a finite time interval At as follows: 


( T Hg,K^ 


new 


~ (Tjjg- g) + 

Mg,K old 


dT 


Hg,K 

dt 


(At) 


where 


dT Hg,K_ v Hg 


2U 7 r TTT 

"dt AX (l ’ H g J K+1 ' ^g.K* + ~ 2 ATaK 


Jp Hg (Cp) Hf 


when (W„ ) = 0, 

Mg liquid 


, 2u z r HT (At)(AT AK ) 
( Hg ’ K new ” ( Hg ’ K old ~ 


J P Hg ( c P) Hg 


In the digital program: 

(1) q k = 2U ? r HT (At)(AT AK ) 
r p Hg^ Cp ^Hg 
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And for the lump containing the interface 


(2) Qj^ - RI x Q k 


ffcond-^.KX 1 - - 


U. 


^ T Hg, K+l T cond 

The temperature at each node K of the NaK and tubes can now be obtained 

5 [a (Cp) NaK p NaK T NaK, k 3 _ 9 C°NaK v NaK p ^ NaK T NaK , k] dA 


at 


ax 


+ U(AT k ) 


HT 


dV. 


NaK 


where a is a factor to add the effect of the tubes, that is, (Cp) e ff ec y ve /(Cp) Na g = a. 
Again assuming P NaK > (Cp) Na K> and are constant and approximating over a 

finite length AX which is equal to the length of one lump, the differential equation 
above can be simplified as 

dT NaK,K _ p NaK v NaK^ Cp ^NaK^ T NaK,K " T NaK,K+P ( 2Ur HT 7r(AX)(AT AK ) 

dt ap NaK^ Cp ^NaK °' W NaK 

40 

where K = 1, 2, 3, . . . , 40. And for a finite time interval At the equation becomes 


(T 


NaK, K 


) . (T NaK K ) ♦ (it) 

new waK ’ K old dt 


where 


dT NaK, K _ v NaK 


dt 


a (AX) ^ NaK > K " T NaK,K+P + 


2U r HT ffL ^ r K^ 


aW 


NaK 


In the digital program 


Q K p Hg( Cp W r27rL _ 2Ur HT L(AT K ) 
aW NaK aW NaK 
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APPENDIX B 


SYMBOLS 


a ht 

(Cp) effective 
(Cp) H g 

^ Cp VaK 

(cp) t 

Ah 

I 

P 

p 

cond 


2 2 

heat-transfer area for condensing region, ft ; m 

2 2 

cross-sectional flow area in mercury condenser tubes, ft ; m 

2 2 

heat-transfer area, ft , m 

effective specific heat of third-loop fluid when tubes are lumped with 
third-loop fluid 

specific heat of liquid mercury, Btu/(lbm)(°F); J/(kg)(K) 

specific heat of NaK, Btu/(lbm)(°F); J/(kg)(K) 

specific heat of tubes in condenser, Btu/(lbm)(°F); J/(kg)(K) 

heat of vaporization for small temperature changes, Btu/lbm; J/kg 

node just before condenser liquid-vapor interface (I is in condensing 
region) 

2 

pressure, psi; N/m 

o 

condensing pressure, psi; N/m 

initial pressure (corresponds to a point on mercury vapor saturation 
curve), psi; N/m 


Q^ n heat flow into condenser from turbine, Btu/(lbm)(sec); J/ (kg) (sec) 

ii. 

change in mercury temperature at k ul node due to conduction which would 
result during time interval At if there were no convection heat flow, 

°F; K 

R specific gas constant, ft-lbf/(lbm)(°R); m-N/(kg)(K) 

RI fraction of lump in condensing region 


r mean inside radius of condenser tubes, ft; m 


r HT 

T 

T 

cond 

T Hg,K 

T NaK 


mean heat -transfer radius of condenser tubes, ft; m 

temperature, °F; K 

condensing temperature, °F; K 

temperature of mercury at node K, °F; K 

NaK temperature, °F; K 
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T NaK, av average temperature of NaK for mercury condensing region, °F; K 
T NaK K NaK temperature at k** 1 node, °F; K 

T q initial temperature (corresponds to point on mercury vapor saturation 

curve), °F; K 

T ri temperature at condenser liquid-vapor interface, °F; K 

(AT)ak average temperature difference between NaK and mercury for one lump 
(K to K + 1) in condenser, °F; K 

AT r temperature difference between NaK and mercury at node K, °F; K 

t time, sec 

At small increment of time, sec 

U overall heat-transfer coefficient (U = U c in condensing region and 

U = Uj in liquid region of condenser), Btu/(sec)(°F)(ft 2 ); J/(sec)(K)(m 2 ) 

U c overall heat-transfer coefficient in condensing region, Btu/(sec)(°F)(ft 2 ); 

J/(sec)(K)(m 2 ) 

Uj overall heat-transfer coefficient in mercury liquid region of condenser, 

Btu/(sec) (°F) (ft 2 ); J/(sec) (K) (m 2 ) 

o q 

V c total volume of mercury vapor in condensing region, ft ; m 

q q 

V NaK volume of NaK in condenser, ft ; m 

Vi total volume of mercury in liquid region of condenser, ft 2 ; m 2 

v Hg velocity of liquid mercury out of condenser, ft/sec; m/sec 

q q 

V Na K velocity of NaK through condenser, ft ; m 

Wnuk weight of NaK in condenser, lbm; kg 

W t weight of tubes in the condenser, lbm; kg 

w Hg mercury liquid flow rate out of condenser, lbm/sec; kg/sec 

w Hg , va po r mercury vapor flow rate into condenser, lbm/sec; kg/sec 

w NaK NaK flow rate through condenser, lbm/sec; kg/sec 

AWjjg difference in flow rates into and out of condenser lbm/sec; kg/sec 

X distance or position, ft; m 


a ratio of effective specific heat of NaK and tubes to specific heat of NaK 

q q 

^Hg specific volume of mercury, ft /lbm; m /kg 
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P Hg 

p NaK. 


Q O 

specific gravity of liquid mercury, lbm/ft ; kg/m 
specific gravity of NaK, lbm/ft^; kg/m^ 
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